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Abstract 

Given a lattice L, a full dimensional polytope P is called a Delaunay polytope 
if the set of its vertices is S D L with S being an empty sphere of the lattice. 



Extending our previous work |DD01] on the hypermetric cone HYPj, we classify 
the six-dimensional Delaunay polytopes according to their combinatorial type. 
The list of 6241 combinatorial types is obtained by a study of the set of faces 
of the polyhedral cone HYPj. 

1 Introduction 

A distance vector {dij)o<i^j<n £ with = ("^^) is called an {n + 1) -hypermetric 
if it satisfies the following hypermetric inequalities: 

n 

H{h)d = bibjdij < for any b = (6i)o<i<„ G Z"+^ with ^ 6^ = 1 . (1) 

0<j<jr<n i=0 

The set of distance vectors satisfying (|1|) is called the hypermetric cone and denoted 

by i/FP„+i. 

In fact, HYPn+i is a polyhedral cone (see ||DGL93[ p. 199). Lovasz (see | peLa97 | 



p. 201-205) gave another proof of it and bound max < n!2"(^^) ^ for any vector 
b = (&j)o<j<n-i defining a facet of HYPn- 

There is a one-to-one correspondence between non- degenerate elements of 
HYPn+i and semi-metrics on affine basis of n-dimensional Delaunay polytopes. So, 
the enumeration of combinatorial types of n-dimensional Delaunay polytopes is re- 
duced to the enumeration of non-degenerate faces of HYPn+i under an equivalence 
relation called geometrical equivalence (see Remark |T| and Section^). Moreover, the 
dimension of faces of HYPn+i allows to define the notion of rank of a Delaunay 
polytope (see |paL92||). 
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We have the inclusion CUTn C HYPn, where CUTn (see Section |^ below and 
Chapter 4 of |PeLa97| |) is the cone generated by all cut semi- metrics on n points. 
One has HYPn = CUTn for n < 6; so, the enumeration of combinatorial types 
of Delaunay polytopes of dimension less or equal to five correspond to the study 
of faces of CUTq. This study was done by Fedorov ( |p^e85|| ), Erdahl and Ryshkov 
P^yEr87|| and [[RyEr88|| ) , and Kononenko ( [|Ko99|| ) in dimension three, four and five. 



respectively. 

In the case of dimensional six, the inclusion CUT-j C HYP-j is strict, but the 
description of facets (i.e. faces of rank 20, corresponding to repartitioning polytopes) 
and extreme rays (i.e. faces of rank 1, corresponding to extreme Delaunay polytopes) 
of the cone HYPj is still possible (see [PD01|| ). In Section ^ we present the number 



of combinatorial types of six- dimensional Delaunay polytope for every given rank. In 
Section ^ we present our method for the computation of those combinatorial types 
using the face lattice of HYPj. 

Voronoi |[Vo08|| defines a partition of the cone PSDn of positive semi-definite 
quadratic forms by L-type domains. Two forms of the same L-type domain have 
affinely equivalent Voronoi polytopes. Every vertex of a Voronoi polytope correspond 
to a center of a Delaunay polytope. All Delaunay polytopes form a partition of the 
space M", which is dual to the partition by Voronoi polytopes. It is proved in [pGL93|| 
that the hypermetric cone HYPn is the union of a finite number of L-type domains. 

While Voronoi theory of L-type domains describes the combinatorial structure of 
lattices, the theory of hypermetrics ( |DGL92| and [PeLa97|| ) describes combinatorial 
structure of one Delaunay polytope in a lattice. While the hypermetric cone HYPn 
has the symmetry group Sym{n), a L-type domain, in general, has trivial group. 



2 Delaunay polytopes and hypermetrics 

Here we present the main notions for hypermetrics, which are needed to our study; 
the presentation is slightly simplified by the systematic use of affine bases. For the 
complete theory, with proofs, see |PGL92|| and Chapters 13-16 of |PeLa97|| . 

A family vq, . . . , f„ of + 1 vertices of is called independent if the family 
{vi — fo)i<i<n has linear rank n. Let L C be a ra-dimensional lattice and let 
S = S{c, r) be a sphere in with center c and radius r. Then S is said to be an 
empty sphere in L if the following two conditions hold: 

(i) ||f — c|| > r for all f G L and 

(ii) the set 5* n L contain an independent set of size n + 1. 



The center c of S* is called a hole in ||CS99|| . The polytope D, which is defined as the 



convex hull of the set S D L is called a Delaunay polytope, or (in original terms of 
Voronoi, who introduced them in [ Vo08|| ) L-polytope. 



Definition 1 Let D be a n-dimensional Delaunay polytope with vertex-set V . 
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(i) A family vq, . . . , Vn of vertices of D is called an affine basis if for all v & V 
there exist an unique family (&i)o<i<n ^ Z""*"-^, such that 



n 



biVi = V and hi = 1 . 

(a) The Delaunay polytope D is called basic if it has at least one affine basis. 
The vertices of an affine basis are called basic vertices. 

All known Delaunay polytopes are basic. We prove in Theorem |1| that all six- 
dimensional Delaunay polytopes are basic too. We will always assume that Delaunay 
polytopes are basic. 

For every family A = {vq, ■ ■ ■ , f m} of vertices of a Delaunay polytope P one can 
define a distance function by dj_{i,j) = \\vi — Vj\\^. The function dj( turns out to 
be an hypermetric by the following formula (see |[As82|| and ||DeLa97|| p. 195) : 

m 

bibjdAit,]) = 2(r2 - \\J2b^Vi - cf) < 0. 

0<i,j<m i=0 

On the other hand, Assouad has shown in | As82 | that every d G HYPn+i can be 
expressed as (i_4 with A being a family of vertices of a Delaunay polytope D of 
dimension less or equal to n. 

A ray d G HYPn+i is called non- degenerate if d = d^ with B an affine basis of 
a n-dimensional Delaunay polytope D. For a given ray d G HYPn+i the annulator is 
defined by 

n 

Ann{d) = {be : = 1 and H{b)d = 0}. 

i=0 

We call basic vectors the vectors Cj = (0, . . . , 0, 1, 0, . . . , 0) G Z"'^^ with < i < n. 
Using Proposition below, those vectors Cj are identified with basic vertices Vi. One 
has H{b) = if and only if 6 = for some < i < n. 

Proposition 1 Let D be an n-dimensional Delaunay polytope with vertex-set V ; let 
d = dji, with A = {vq, . . . ,Vn} a family of V. Then there is equivalence between 
following properties: 

(i) the family A is independent and 

(li) det [d{0,i) +d{0,j) - d{i, j)]i<ij<n 7^ 0. 

Also one has equivalence between properties: 

(Hi) A is an affine basis and 

(iv) the following mapping is one-to-one 



b G Ann{d) t-^ ^ biVi G V . 



1=0 
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Proof. The family (fi)o<j<n is independent if and only if the matrix ((fj — vo).{vj — 
VQ))i<ij<n is positive definite. The formula 2(t>j — Vo).{vj — vq) = \\vi — fo|p + \\vj — 
foll^ — \\vi — VjW^, gives the first equivalence. The second equivalence is obvious. 

By above Proposition, the set Ann{d) is finite if d is non-degenerate, since the 
hypermetric vectors b G Ann{d) correspond to lattice points v, which belong to a 
sphere. We will consider only non-degenerate hypermetrics (in the hypermetric cone 
HYPn+i) associated with affine bases of n-dimensional Delaunay polytopes. 

Using above mapping, we identify vertices v of D with hypermetric vectors b. 
For example, for the affine basis B = {vq, . . . ,^3} with vq = (0,0,0), vi = (1,0,0), 
V2 = (0, 1, 0) and = (0, 0, 1) of the 3-cube, one has (1, 1, 0) = f 1 +f2 — fo, (0, 1, 1) = 
V2 + V3 — Vq, (1, 0, 1) = f 1 + f3 — Vq and (1, 1, 1) = f 1 + f 2 + f 3 — 2t>o. Note that, while 
(—2, 1, 1, 1) G Ann{ds), the inequality H{—2, 1, 1, l)d > does not define a facet of 
the cone HYP^, since H{-2, 1, 1, 1) = H{-1, 1, 1, 0) + H{-1, 0, 1, 1) + H{-1, 1, 0, 1). 

The face F{d), associated to an hypermetric d G HYPn+i, is the minimal face, 
containing the vector d. It can also be defined as: 

F{d) = {ee HYPn+i : H{b)e = for all b G Ann{d).} 

The rank of a n-dimensional Delaunay polytope D is defined as dimension of 
F{dis), where B is an affine basis of D (see ||DeLa97|| p. 217); its corankis, by definition, 
("2^^) ~ rank D. The rank of a Delaunay polytope D is equal to the topological 
dimension of the set of affine bijections T of M" (up to translations and orthogonal 
transformations), for which T{D) is again a Delaunay polytope (see [peLaQ?]] p. 225). 



If D is an n-simplex, then its unique affine basis defines an hypermetric d, which 
has Ann{d) = {eo, . . . , e„} and F{d) = HYPn+i- If -D is a n-dimensional Delaunay 
polytope with n + 2 vertices Vq, . . . , Vn+i (such polytope is called repartitioning poly- 
tope) and affine basis B = {vq, . . . , t>„}, then, by writing Vn+i = Yl^=o ^i'^i^ obtains 
Fide) = {d e HYPn+i,H{b)d = 0}. So, F{d) is a facet of HYP^+i. On the other 
hand, if G HYPn^i with F{d) being a facet, then Ann{d) = {cq, . . . , e„, 6} and 
d = ds with B being an affine basis of a repartitioning polytope. 

Define two Delaunay polytopes to be affinely equivalent, if there is an affine 
bijective mapping transforming one into another; the equivalence classes by this rela- 
tion are called combinatorial type. By Theorem 15.2.1 in [ DeLa97|| p. 222, two vectors 



d, d' of HYPn+i, such that F{d) = F{d'), correspond to affinely equivalent Delaunay 
polytopes. Moreover, if two Delaunay polytopes D, D' are affinely equivalent, then 
an affine basis B of D induces an affine basis B' of D', such that F{ds) = F^dsi). 

A face of HYPn+i is called non- degenerate, if one of its interior vectors d is non- 
degenerate (by Theorem 15.2.1 of | DeLa97|| , this is equivalent to all interior vectors 



being non-degenerate). One can consider only non-degenerate faces, since we already 
assumed, that all Delaunay polytopes are basic. Given a face F of HYPn+i, we define 

Ann{F) = {6 G Z"+i : ^ 6^ = 1 and H{b)d = for all d G F}. 

i 

We will call {6°,..., 6"} C Ann{F) an affine basis of the combinatorial type if 
det {W , . . . yb"^} = ±1. Any affine basis of Ann{F) corresponds to an affine basis 
of a Delaunay polytope D, if D has this combinatorial type. 
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Definition 2 Two faces F, F' of HYPn+i are said to be geometrically equivalent if 
there exist an affine basis B — {6°, . . . , 6"} of Ann{F), such that the mapping 

f 4)F',F ■■ Ann{F') Ann{F) 

\ ' b ^ 6o&° + ...M" 

is bijective. 

If TT is a permutation of {0,...,n} and F a face of HYPn+i, then 
{e7r(o), • • • ! e7r(n)} IS an afiinc basis of Ann{F)] so, two faces F and F', which are 
equivalent, by a symmetry of HYPn+i, are also geometrically equivalent. The reverse 
is not true, in general. 

Remark 1 Let D be a n- dimensional repartitioning polytope with an affine ba- 
sis B — {t'o,...,t'n}; we write = ^"=o^«'^i- family B' — 
{vq, . . . Vi-i, Vn+i, fi+i, . . . , Vn) is an affine basis if and only if \bi\ — 1. 

If bi = -1, then F{dB') = Fid^) = {d E HYPn+i,H{b)d = 0}, while 
if bi = 1, then F^ds) ^ F{dB') ^ {d e HYPn+uH{b')d = 0} with b' = 
{—bo, . . . , —bi-i, bi, —bi-i, . . . , —bn)- 

So, the study of combinatorial types of n-dimensional Delaunay polytopes is 
reduced to the study of non- degenerate faces of the hypermetric cone HYP^j^i under 
the geometrical equivalence. 

Consider, as an example, the two-dimensional case. We have HYP^ — CUTs 
and CUTs has three facets, which correspond to H{— 1,1,1) and its permutations. 
One can check that a triangle T = {vq,vi,V2} satisfies dr = {doi,dQ2,di2) € HYP^ 
if and only if it is an obtuse triangle. The vector dq- is non-degenerate if and only if 
three vertices of T are not aligned. Moreover, dr is incident to an hypermetric facet, 
say, H[— 1,1,1), if and only if the vertex vq has angle |, in which case the Delaunay 
polytope has four vertices: Vq, Vi, V2 and Vi-\-V2 — Vo. So, there are two combinatorial 
types of Delaunay polytopes in dimension two: obtuse triangles and rectangles. 

Take a face F of the hypermetric cone HYPn+i- To every b = (6i)o<i<n G 
Ann{F), we associate a vertex {bi, . . . ,bn) G 'Z/^', let the set of such vertices being 
denoted by V. Every distance vector d E F correspond to a Gram matrix G. The set 
V is the vertex-set of a Delaunay polytope, which is described by the scalar product 
defined by G. So, we can encode, in our computations, the combinatorial type of 
Delaunay polytope by the set Ann{F). 

Proposition 2 If D is a Delaunay polytope andB, B' are two affine bases ofD, then 
the faces F{dB), F{dB') are equivalent up to a linear mapping. This linear mapping 
preserve the non- degeneracy. 

Proof. If {vq, . . . ,Vn) and {v'q, . . . ,v'^) are two affine bases of D, then one can express 
in terms of {vj)o<i<n as follows 

n n 

v'^ — OiijVj with aij — 1 . 
j=o j=0 
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One can express | |f — v'^^ | P in terms of | |f ~ "^j^ I This induces a hnear mapping 
from F^dis) to F{ds'); expressing Vj in terms of (fj')o<i<n, one get the reverse mapping 
(^~^ and so, the hnear equivalence. 

If d E F{dis), then d is non-degenerate if and only if (j){d) is non-degenerate. 

Definition 3 Let B be an ajfine basis of a Delaunay polytope D; then 
(i) If rank D = 1, then D is called extreme. 

(a) If all sub-faces of F^d^) are degenerate, then D is called maximal. 

Above definition of maximality is independent of the choice of affine basis since, 
by above Proposition, the linear equivalence between two faces preserves the non- 
degeneracy. Obviously, any extreme Delaunay polytope is maximal. We present in 
Corollary ^ the list of all maximal six-dimensional Delaunay polytopes. 

Let 7n = {0, 1}" be the vertex-set of the Delaunay polytope of the lattice Z" 
and let 

n 

/i7„ = {x G 7„ : Xi is even} 

i=l 

be the vertex-set of the Delaunay polytope with center c = (1/2, . . . , 1/2) (this poly- 
tope is called half n- cube) of the root lattice 



F>n = {x G Z" : "^^^Xi is even} . 



1=1 

Proposition 3 The n-cube 7„ and the half n-cube h'jn have rank n and are maximal 
Delaunay polytopes. 

Proof. One can define the following distance functions on 7^ and /i7„: 

c^i : 7n X 7n 

(x, x') (— > (Xj — x-)^- 

If d is the distance function of the n-dimensional Delaunay polytope 7„ or /i7„, then 
this distance function is expressed as 



d = Xidi with Xi > 0, 



i=l 



which proves that the rank of 7^ and of /i7„ is n. Now, if a face F is included in -F(7n) 
or F{h'jn), then one of the coefficients A, becomes zero and the face F corresponds 
to 7„_i or hjn-i, which is of lower dimension n — 1. So, the Delaunay polytopes 7„ 
and h'jn are maximal. 
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3 The case of dimension six 



Call cut cone and denote by CUTn+i the cone generated by all cuts 6s G (where 
5" is a subset of {0, ... , n}), defined by 



(5, 



S)ij 



1 if \Sn{i,j}\ = 1 and (5, 



S)ij 



0, otherwise. 



Clearly, H{b)6s = b{S){l—b{S)) with b{S) = J2aes ^^^^ proves that all hypermetric 
inequalities are valid on CUTn+i. So, C[/T„+i C HYPn+i- Moreover, a cut Ss is 
incident to the face, defined by H{b), if and only if b{S) = or 1. 

The list of 3773 facets of HYP^ was found by Baranovskii ||Ba99|| using the 
method described in [[Ba70|| , i.e. he found, by hand, that for all other hypermetric 
vectors 6, one can express H{b) as a sum of terms H{b') with b' belonging to his list of 
3773 elements. While this result was announced in |[Ba99|| , the detailed computations 
were not published. In ||DD01|| , another method was proposed: if the Baranovskii's 
list was not complete, then, in our computation |pD01|| of the extreme rays of HYP-j, 
we should find some extreme rays, which are not hypermetric. But this was not the 
case; so, the list is complete. 

The list of representents of 14 orbits of facets is given below: 



b' = 


(1,1,-1,0,0,0,0) 




62 = (1,1,1,-1,-1,0,0) 


b' = 


(1,1,1,1,-1,-2,0), b" 


= (2,1,1,-1,-1,-1,0) 


b' = 


(1,1,1,1,-1,-1,-1) 




fe6 = (2,1,1,1,-1,-1,-2) 


b' = 


(2,2,1,-1,-1,-1,-1), 




= (1,1,1,1,1,-2,-2) 


b^ = 


(3,1,1,-1,-1,-1,-1), 




= (1,1,1,1,1,-1,-3) 


6" = 


= (2,2,1,1,-1,-1,-3), 


b^'' = 


= (3,1,1,1,-1,-2,-2) 


fei^ = 


= (3,2,1,-1,-1,-1,-2), 


b'' 


= (2,1,1,1,1,-2,-3) 



It gives the total of 3773 inequalities. The first ten orbits are the orbits of 
hypermetric facets of the cut cone CUT-i] first four of them come as 0-extension 
of facets of the cone HYPq = CUTq (see |PeLa97|| , Chapter 7). Last four orbits 



consist of some 19-dimensional simplex-faces of CUTj, becoming 20-dimensional, i.e. 
simplex-facets in HYPj. 

Using Remark |l], we obtain that the list of 14 orbits of facets fall into nine 
equivalence classes b\ b\ 6^ 6^ {&^&^}, {b^\b^^}, {b^^,b^^}. So, 

there are nine combinatorial types of six-dimensional Delaunay polytopes of rank 20 
(i.e. repartitioning polytopes). 

In |PeLa97|| p. 229 another notion, called switching by root of b, is defined: if 
A C {0, . . . ,n} and b{A) = 0, then define 6^ by bf = -bi ii i e A and bf = h if 
i ^ A. It is proved that the switching by root of a facet- defining vector of HYPn+i is 
again a facet- defining vector. Three following vectors define facets of HYPj, which 
are switching by root equivalent: 

b' = (2, 2, 1, -1, -1, -1, -1), b' = (2, -2, 1, 1, 1, -1, -1), b' = (-2, -2, 1, 1, 1, 1, 1). 

On the other hand, b^ and b^ are geometrically equivalent by Remark |l], while b^ and 
b^ are not geometrically equivalent. 
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Remind that Eq, E7, Eg are root lattices defined by 



Eq = {x e El : xi + X2 = 0} , Et = {x e Es : Xi + ■ ■ ■ + Xs = 0}, 
Es = {x eM.^ : X eZ^uU+ and ^. Xi e 2Z} . 



Unique type of Delaunay polytope of Eq is called Schlafli polytope (see ||Cox63|| ) 



and denoted by Sch. Its skeleton graph is 27-vertex (strongly regular) graph, called 
the Schlafli graph, whose symmetry group has size 51840 and is isomorphic to 
the group of isometry, preserving the Schlafli polytope. This group is denoted by 
Aut{Sch). 

In |PGL92|| were found 26 orbits of non-cut extreme rays of HYPj by classifying 



the afiine bases of the Schlafli polytope of the root lattice Eq. For every non-cut 
extreme ray (K+f) of HYP^, there exist a facet-inducing inequality f{x) > of 
CUTf, which is non-hypermetric, so that f{v) < 0. This property establish a bijection 
between the 26 orbits of non-hypermetric facets of the cut cone CUT^ (see ||Gr9CI|| ) 
and the 26 orbits of non-cut extreme rays of HYP^ and proves that HYPj has 29 
orbits of extreme rays: three orbits of non-zero cuts and 26 orbits coming from Sch 
jDDOTl ). 



see 



Proposition 4 Let B be an affine basis of Schlafli polytope; then 

(i) The distance vector is incident to 20 hypermetric faces of HYP-j, which 

are all facets of HYPj. 

(a) If F is a face of HYPj, containing the vector d^, then it is non- degenerate 

and \Ann{F)\ = 7 + corank{F). 

Proof. The Schlafli polytope is six-dimensional and has 27 vertices. So, for every 
affine basis B, the vector d^ satisfies H{h)ds = for 27 different b G Annld^)] we 
write Ann{dQ) — {eo, ■ ■ ■ , e^} = {b^, . . . , 6^"}. On the other hand, it is known ( |pGL92 | 



and |PeLa97|| p. 239), that the Schlafli polytope has rank 1. So, the rank of the matrix 
[/7(6i),...,i7(620)] must be {^f) - 1 = 20. So, the family [H{b%<i<20 is linearly 
independent. If one of H{U) is not a facet of HYPj, then it can be expressed in terms 
of [HiVyij-Li, this contradicts to linear independence and so, (i) holds. 

The extreme ray ds is non-degenerate and is a sub-face of F; so, F is also 
non-degenerate. Every hypermetric face, containing F, contains ds'-, so, one has 
Ann{F) — {cq, . . . , eg} = {b'^ , . . . , b'^} C {6^, . . . , 6^°}. The linear independence of 
the family [-&(&*)] i<j<20 implies that k = corank{F) and so, (ii) holds. 

Above Proposition is not true for the 56- vertex Gosset polytope (| |Gox63| |): the 
Gosset polytope has 374 orbits of affine bases. Each extreme ray, corresponding to an 
affine basis {vq, . . . , ^7} of the Gosset polytope, is incident to 48 (= 56—8) hypermetric 
faces of HYPg,. But amongst these 48 face-defining inequalities, the number of facets 
varies from 27 (= (^^^) — 1) to 41. See ||DeLa97| j , p. 230 for general lower bounds (on 
the number of vertices of a Delaunay polytope) as a function of its rank. 

Theorem 1 All six- dimensional Delaunay polytopes are basic. 
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Proof. The simplex is a basic polytope, since there are 7 vertices and they form an 
affine basis. 

Assume that D is a non-simphcial Delaunay polytope of a lattice L generated 
by the vectors wi, . . . , wq. Denote by V the volume of the simplex formed by the 
vectors 0, wi, . . . , w^. Take a family of 7 independent vertices A = {vq, . . . , vq} in 
the vertex-set of D and denote the volume of the corresponding simplex by V. 

In [ [RyBa98|] , it was proved that the relative volume = ^ is 1, 2 or 3. If A; = 1, 
then A is an affine basis and we are done. Assume now that k > 1, i.e. that A is 
not an affine basis. Then, there exists a vertex v of D, which is written uniquely as 
V = Yl^=o ^i'^i with b being fractional and Yl^=o = ^■ 

The distance vector c?^ satisfies H{h)dj\^ = with b = {bo, ...,66) being a frac- 
tional hypermetric vector. But one can express H{b) as XiH{b^) with > 
and fe' being a permutation of one of the following vectors (see |PR,yBa98| ) : 



Case k = 2 


CcisG k — 3 


i(-l,-l,l,l,l,l,0) 
i(-l, -1,-1, 1,1, 1,2) 
|(-2, -1,-1, 1,1, 1,3) 
i(-2,-l,l,l,l,l,l) 
i(-l, -1,-1, -1,1,2,3) 
i(-3,-l,l,l,l,l,2) 
(-1,1,1,0,0,0,0) 


i(-l, -1,-1, 1,1,2,2) 
i(-l, -1,-1, 1,1, 1,1) 
|(-2,-l,l,l,l,l,2) 



Since A/ > 0, one has H{b'')d = 0, i.e. the vectors w'l = X^Lo ^i'^i ^ith I < I < N 
are vertices of the Delaunay polytope D. 

Since b is fractional, at least one of 6' is fractional, say, 6'°. But all fractional 
hypermetric vectors of above Table have one coordinate with absolute value equal to 
|, say, = ^. So, the family {t>o, . . . , f^'j^, fi+i, . . . , fe} defines a simplex of 
relative volume = 1, i.e. it is an affine basis. 



Theorem 2 The 6241 combinatorial types of Delaunay polytopes are partitioned by 
rank in the following way: 
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rank 


Nr. in MY 


Nr. m Uu I7 


21 


1 (simplex) 





20 


9(repart.) 


1 


19 


30 


2 


18 


95 


8 


n 


233 


28 


16 


500 


95 


15 


814 


241 


U 


1092 


434 


13 


1145 


527 


12 


984 


481 


11 


686 


325 


10 


417 


183 


9 


218 


83 


Q 

Q 


1 Uo 


00 


7 


52 


13 


6 


21 


3 


5 


8 





4 


4 





3 


2 





2 


1 





1 


l(Schlafli) 






Proof. The proof is purely computational and the method is described in next 
Section. 



Corollary 1 All maximal six- dimensional Delaunay polytopes are: Schlafli polytope, 
6-cube, half 6-cube and direct product of half 5- cube with 1-cube. 

Proof. This result follows directly from the computation of above Theorem. 



4 Computational methods 

Our computation of combinatorial types of six-dimensional Delaunay polytopes used 
the face-lattice of HYPj; combinatorial types of Delaunay polytopes of corank i + 1 
were found from combinatorial types of Delaunay polytopes of corank i. We start 
from the list of combinatorial types of corank 1, i.e. the nine repartitioning polytopes. 
The plan of our computation was as follows: 

(i) Take the list of combinatorial types of Delaunay polytopes of corank i. 

(ii) For each of them, find all sub-faces, using our knowledge of facets and extreme 
rays of HYPj; we obtain faces of corank i + 1. 

(iii) For every face F of corank i -|- 1, find extreme rays (M+/j)i<j<Ar, contained in it, 
and define d = X^ili /«• The distance vector d is in the interior of F; so, using 
Proposition |I|, one can test if this ray is non-degenerate or not and this tells us 
if the face is non-degenerate or not. 

(iv) Find the classes of geometrical equivalence amongst the non-degenerate faces 
and so, the list of combinatorial types of corank i + 1. 
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Above procedure finds all combinatorial types of Delaunay polytopes from 
corank 1 (i.e. repartitioning polytopes) till corank 20 (i.e. the Schlafii polytope). 
But in order to describe completely our method, we need to precise how we find the 
classes under geometrical equivalence. 

The first algorithm for the problem of geometrical equivalence is the following: 
given two faces F and F', find all affine bases of Ann{F) until one finds an equivalence 
4>F,F' (see Definition |]). This method works for corank 1 or 2 and was used by 
Kononenko in the five- dimensional case. But the number of affine bases becomes too 
important to be workable in corank 3. 

So, one needs another, more efficient method. The first idea is to split the set of 
faces of HYP-j into two classes: faces, which contain (one or more) Schlafii extreme 
rays, and those, which are generated only by cuts. 

Definition 4 Let F he a face of HYP^, which contains a Schlafii extreme ray 
corresponding to an affine basis B = {vq, . . . ,vq}. Then, every b G Ann{F) defines a 
vertex v = X]i=o^«^* of Sch. All such vertices are denoted by S{F,dB)- 

So, every combinatorial type of faces, which contains at least one Schlafii ex- 
treme ray, can be interpreted as a set of vertices of Sch. 

Proposition 5 Let F and F' be two faces, which contain at least one Schlafii extreme 
ray; then it holds: 

(i) If F and F' are geometrically equivalent, then for every Schlafii extreme ray 
ds in F, there exist an affine basis B' of Sch, such that S{F,dB) = S{F',ds/). 

(a) If F and F' contain a Schlafii extreme ray d^ and d^i, such that S{F, d^) is 
identical to S{F', d^/) up to an element of Aut{Sch) , then F and F' are geometrically 
equivalent. 

Proof. For every distance vector dis G F, one gets an identification of Ann{F) with 
S{F, djs). If F' is geometrically equivalent to F, then the vectors Cq, . . . , eg in Ann{F') 
are identified with vertices Vq, . . . ,Vq in S{F, d^). These vertices form an affine basis 
B' and one has S{F' , d^i) = S{F, dts); so, (i) holds. 

Since S^F.djs) is isomorphic to S{F',dj3') by an element of Aut{Sch), one can, 
without loss of generality, assume that they are identical. Now, B is identified with 
vertices vq, . . . ,vq of S{F, d^) and B' is identified with vertices v'q, . . . ,Vq of S{F, d^). 
Then, the expression of v'^ in terms of Vi determine an affine basis of Ann{F), for 
which the mapping is well-defined and bijective. 

Above Proposition express the geometrical equivalence in terms of the existence 
of an element of Aut{Sch) mapping a set of vertices into another set of vertices. 
Those sets of vertices are identified with corresponding sets of vertices of the Schlafii 
graph (the Schlafii graph and the Schlafii polytope have the same symmetry group 
Aut{Sch)). So, the problem is expressed in graph-theoretic terms and can be solved, 
using, for example, the nauty program ( ||MK| | ) . Therefore, one can build the geomet- 
rical equivalence classes. 
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Now, we extend above method to the case of faces generated by cuts. Let F be 
a face, generated by cuts {S{Si)}i<i<iy; we first need to find Ann{F), i.e. all vectors 
b e Z"+i with Y.Lo bi = 1, having H{b)6{Si) = 0. 

Those equations can be rewritten as '^^eSi = Xi with Xj = or 1, i.e. a linear 
system in b. This linear system has rank n + 1, because of Proposition below; so, one 
can find the set Ann{F) for every face, generated by cuts. 

Proposition 6 Let T he a face of HYPn+i generated by cuts (5(5'.t))i<i<Ar. Then the 
following properties are equivalent: 

(i) the face F is non- degenerate and 

(a) the linear system, formed by the equations Ylxes^^ ~ ^ ^'^^ Y27=o^i ~ ^' 
has solution set {0}. 

Proof. If the face F is degenerate, then there exist a vertex v, which can be expressed 
in two different forms v = X]r=o ^i'^i ~ Sr=o ^i"^*- lifter denoting ai = bi — 6-, one 
gets V = X]r=o(^* + kai)vi with k E Z and b + ka belongs to Ann{F). Therefore, one 
gets (6 + ka){Si) = or 1, and Yl^=o + ken = 1 for all e Z. This is possible only 
if Exe5, = and ^^^^ = 0. 

Let the solution set be non-zero, i.e. suppose that one can find an integer-valued 
non- zero solution a. This implies that the vectors CQ + ka belong to Ann{F) for every 
A; G Z. So, Ann{F) is infinite and F is degenerate. 

If we write Ann{F) = {b^, . . . , b^}, then every cut 6{Si) G F with 1 < i < 
defines an Euclidean semi-metric on the set {eo,...,e„}; this semi-metric can be 
uniquely extended to Ann{F) by 6i{b,b') = \b{Si) — b'{Si)\. 

So, to every face, generated by cuts, one can associate a set of semi-metrics on 
Ann{F), which are, in fact, cut semi-metrics. 

A combinatorial type of Delaunay polytope encodes all possible embeddings of 
Ann{F) into the vertex-set of a Delaunay polytope of a lattice. These embeddings 
are completely described by the distance vector dy on their vertices. This distance 
vector is expressed as J2iLi ^i^i with Aj > 0. 

Therefore, the combinatorial type of a face, generated by cuts, corresponds to 
the description of all semi- metrics on this set. This information can be expressed in 
graph-theoretic terms; so, we can test if two faces, generated by cuts, are isomorphic, 
using the nauty program ( ||1V1K|| ) . So, again one can build the geometrical equivalence 



classes and our method is completely described. 
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